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A many-body wavefuction is postulated, which is sufSciently general to describe 
superconducting pair-correlations, and/or spin-correlations, which can occur either 
as long-range order or as finite-range correlations. The proposed wave-function 
appears to summarize some of the more relevant aspects of the rich phase-diagram 
of the high- Tc cuprates. Some of the states represented by this wavefunction are 
reviewed: For superconductivity in the background of robust anti-fcrromagnetism, 
the Cooper-pairs are shown to be a superposition of spinquantum numbers S=0 and 
S=l. If the anti-ferromagnetism is weak, a continuous super-symmetric rotation 
is identified connecting s-wave superconductivity to anti-ferromagnetism. 



1 Introduction 

The large variety of phenomena due to correlations between electrons in solids 
is one of the major challenges of contemporary physics. Over the past 80 
years this subject has formed the arena of many new developments in theo- 
retical physics. The most common approach to the many-body problem is, 
to define in the first step a model defined as a Hamiltonian, Lagrangian, or 
action, describing the most important interactions of the system one tries to 
understand. In the second step one tries to solve the corresponding equations 
of motion. Apart from a few special cases, no exact solution is known of the 
Schrodinger equation of a large number of interacting particles. To circum- 
vent this problem two different approaches are most frequently used, often in 
combination with each other: One is to device an analytical approximation 
scheme, which ultimately may provide an approximate solution. The other is 
to use computational techniques for a finite size cluster. 

An alternative approach, which has sometimes been quite successful, is to 
start with designing a variational many-body ground-state wave-function with 
the desired characteristics of (for example) the pair-correlation function, spin- 
polarization, etcetera. One can than subsequently try to find a model Hamil- 
tonian which will produce (i) the aforementioned many-body wave-function as 
it's ground-state, and (ii) the low energy excitations which are important for 
determining the thermal properties and the various spectral functions all of 
which can be measured experimentally. 
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The variational wave-function approach has been used with great success 

for superconductivity^, the fractional quantum Hall effect?, and the resonat- 
ing valence boncf Other striking examples of this 'bottom-up' approach 
are the Affleck-Marsten flux-phas^'^'^°'^^'^^'^^, or d-density wave^^, the mixed 
anti- ferromagnetic/superconducting state^^, and the Gossamer superconduct- 
ing wavefunction'^^'^^. 

Here we use the many-body wave-function approach to explore the possi- 
bility of (anti)-ferromagnetism, superconductivity, and local pairs within the 
same phase-diagram. Our ansatz for the many-body wavefunction is 

|*) = Cexp{5]5^/(fc,G) 

G k 

where C is a normalization constant. The function /(fc, G) is a distribution 
function of the center of mass momentum (G) and the relative momentum (2fc) 
of pairs of electrons. Below we will see, that /(fc, G) represents several types of 
different correlations within the interacting electronic system. Moreover Eq. 1 
is sufficiently general to capture, besides less familiar states of matter, several 
well known states of matter, such as superconductivity, anti-ferromagnetism, 
resonating valence bond states, and the non-interacting electron gas. It there- 
for may be a useful starting point for the rich phase-diagram of the high-T^ 
cuprates, and of the heavy- fermion superconductors. In the subsequent dis- 
cussion we will frequently refer to the projection of Eq.l on the subspace with 
27V^ electrons 

|*2iv) = /(^' G')4ct_,^c J^|0) (2) 

G k 

This expression is a generalization of the many-body wavefunction describing 
the simultaneous occurrence of superconductivity and anti-ferromagnetism^^ 
to arbitrary distributions, f{k, G), of the center of mass momentum G. 

2 Superconductivity 

Let us first consider the special case, that the distribution function f{k,G) is 
a Dirac J-function of the center-of-mass momentum: f{k, G) = 6{G, 0)'^/j. 

14-2^) =CP2Jvexp{4}|0) (3) 

where P2JV is the projection on states with 27V electrons, and 

4 = E'^4tcU (4) 

k 
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creates a pair of electrons with zero center-of-mass momentum. The exponen- 
tial function of appearing in Eq.3, can be expanded in a Taylor's series. 
The term containing 2N electrons 

|*2A.)={4}^|0) (5) 

is similar to a Bose-Einstein condensate of N composite bosons, each con- 
sisting of a pair of electrons. In a superconductor this state is realized if we 
cool an isolated lump of material with known number of electrons below the 
phase transition. Bardeen, Cooper and Schrieffer considered a phase-coherent 
superposition of states containing different numbers of pairs: 

|vl>^)=C^e^<^^|vE'2^) (6) 

JV 

which can be shown to be equivalent to the state 

\'^4>)=Jl{uk + e">'vkcl^cl^^)\Q) (7) 

k 

where Vfc = Vk/uk, and \vk\'^ + jitfep = 1. In the limit case, where jufep 1 
for |fc| < kp and \vk\'^ — > for > kp, this product is just the ground-state 
of the non-interacting electron gas 

i*o)=nc4jo) (8) 

where the occupation number 9^ is or 1 depending on whether the state is 
above or below the Fermi level respectively. 



3 Anti-ferromagnetism 

Another well-documented limit of Eq. 2 is presented by the anti-ferromagnetic 
state. To be specific, we consider a square lattice in two space-dimensions, 
with site-alternating up- and down polarization. The anti-ferromagnetic Bragg 
vector is Q=(7r, tt) in this case, and the operators which create an electron on 
the two sublattices are 



kcr — ^ y-ka H- 
"fecr — ^ \^ka ^ ^k- 



(9) 



-'k+Qa 
-^k+Qa 

respectively. The single particle operators of the anti-ferromagnetic state are 
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where the coefficients satisfy /Ufc = /i_fc, Uk = v-k, l^k+Q = i^k, t^k+Q = Hk, 
and l/ifep + \i^k\'^ = 1- The energies are determined by the details of the 
energy-momentum dispersion of the c|.^-operators and by the interactions giv- 
ing rise to the anti-ferromagnetism. For simpHcity we consider here the case 
where only one of the two sub-bands is occupied, which assumes that the 
anti-ferromagnetic splitting is very large. The ground-state wave-function is 
obtained by partially or fully occupying the lowest band 

i*AF)=nc4T4iio> (11) 



The states with momentum k and —k are degenerate, so we may replace al.-^f3l^ 

with Q!fci-/3!Lfcj^ in the above product. Because each electron-state can only be 
created once, this is equivalent to 

\^AF) = {2-'/'J2^^ayU,}^\0) (12) 

k 

The derivation of Eq. 12 employs the property that aj^^ and /3^^ are orthogonal 
to each other, consequently all two-particle operators under the summation 
commute with each other. In the first (reduced) Brillouin-zone of the anti- 
ferromagnetic state the occupation function 6^ is either 1 or 0. In Eq. 12 
and in later similar expressions the summation over k refers to the paramag- 
netic (extended) Brillouin-zone. The factor 2^^/^ compensates for the double 
counting. However, aj,^ = —(^I-^qu a-nd pj.^ ~ /3fc_|_Qg.- To avoid the complete 
cancellation of terms originating from the first and second reduced Brillouin- 
zone, it is therefor important to define 0^_^q = —0^ for k in the second reduced 
Brillouin-zone. 

To see the relation with Eq. 2 we substitute for and the origi- 

nal c^^-operators in the paramagnetic Brillouin-zone using Eq. 10. We then 
decompose the summation over k in separate terms corresponding to different 
quantum numbers for the spin and the center-of-mass momentum 



= Sfe /o(^)4t'^-/s1 + J^k'id'^ik') +S'°(fc))4'-|-Q/2-r^-fe'-|-Q/2J. 



(13) 



where 

fS{k)= 0f(2^fcz.fe) 

€+Q/2{^^l'+Q/2-''l'+Q/2) (14) 
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This concludes the proof, that Eq. 2 is sufficiently general to also con- 
tain the anti-fcrromagnctic state as one of the possibiUties. The distribution- 
function of the center-of-mass momentum, f{k,G), is a superposition of two 
^-functions: The first at G = is an even function of k and therefor cre- 
ates a pair with spin quantum number zero. The second 5-function at the 
anti-ferromagnetic Bragg-vector G = (tt, tt) contains two contributions with 
different spin quantum-numbers. Since 6-k'+Q/2 = ^fc'-Q/2 = —^k'+Q/2, and 
^-fe'+Q/2 ^ ^■'^k'-Q/2 ^ ^'k'+Q/2- ^^'^ functions cf{k') and g°{k') are even and 
odd functions of k' respectively. Hence the first of the two terms with center- 
of-mass momentum Q corresponds to a S=0 state. The other term, because it 
is an odd function of k' , corresponds to the ms = member of the {S = 1) 
spin-triplet. 

A well-known property of the anti-ferromagnetic state is the broken SU(2) 

symmetry. This symmetry breaking is illustrated by the fact, that Eq. 13 cor- 
responds to a superposition of pair-states with different spin quantum numbers 
{i.e. S=0 and S=l). Eq. 13 is also a superposition of pair-states with differ- 
ent center-of-mass quantum numbers, and this illustrates the fact that the 
anti-ferromagnetic state also breaks the discrete translational invariance of the 
lattice. 

The spin-polarization depends on the value of /i^ = 1 — i^^. For example 
Mfe = ^ft = l/V^ corresponds to full spin-polarization, where the electrons 
move either completely in sublattice a with spin |, or in sublattice b with spin 
J,. From Eq. 14 we sec that the fully polarized anti- ferromagnetic state can be 
regarded as a simultaneous condensation in a singlet state with momentum G = 
and a triplet with momentum G = (tt, tt), both having the same amplitude. 

4 Superconductivity and anti-ferromagnetism 

The results of this section and the previous one can be summarized by a dia- 
gram displaying the center-of-mass distribution function f{k, G) averaged over 
the relative coordinate k, F{G) = |/(/e,G)p. In Fig. 1 this function is 
displayed for the superconducting state and the anti-ferromagnetic state. We 
point out two important aspects of the BCS-wavef unction of the superconduct- 
ing state: 

1) The internal structure of the Cooper-pairs {i.e. the pairing symmetry 
and the degree of localization of the pairs) is completely determined by Vk 
and Uk in Eqs. 6 and 7: The function tpk = Wfe/wfe is the Fourier-transform of 
the wavefunction describing the relative coordinate of the electrons forming a 
Cooper-pair^^. 

2) All pairs are condensed in a state with zero center-of-mass momentum 
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(0,0) Q (7t.7c) (0,0) Q (n,n) 

Figure 1: Distribution of the ccntcr-of-mass momentum for a) a. BCS-superconductor b) an 
anti-fcrromagnct. The positive (negative) values correspond to the singlet (triplet) channels. 
The vertical bars indicate (5-functions, the area of which is represented by the length of the 
bar. 

G = 0, regardless of the details of the internal structure of the pairs. 

Although the long-range anti-ferromagnetic state is characterized by two 
^-functions for the center-of-mass momentum G, Eq.l3 does not correspond to 
a true superconducting state, because the occupation function 9^ has a sharp 
cut-off at the Fermi-momentum. 

4-1 Superconductivity in the background of anti-ferromagnetism 

We first consider the situation, where the anti-ferromagnetism has been sta- 
bilized on a high energy scale, and superconductivity is a relatively weak 
phenomenon. Under those conditions it is a good approximation to assume 
that the onset of superconductivity does not affect the anti-ferromagnetic or- 
der parameter. The simtiltancous occurrence of superconductivity and anti- 
ferromagnetism requires that 6^ is replaced with the smooth function Vk/uk, as 
usual in a BCS-type superconductor. In principle this smooth function may be 
different for the three different terms with differing quantum numbers, and the 
corresponding many-body wavefunction was described by Chen et all^ , but in 
the fully polarized anti-ferromagnet we expect a single distribution function, so 
that the supcrconducting/anti-ferromagnetic wavefunction can be transformed 
to up- and down-spin sublattice-operators using Eq. 9 

\^AF) = \['^ai,hi,\0) (15) 
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Interestingly in this example the Cooper-pairs are a superposition of an S=0 

and an S=l spin-state with equal amplitudes for the S=0 and S=l contri- 
butions. In other words, the superconducting state of a fully polarized anti- 
ferromagnet is a condensate with a broken SU(2) symmetry, with Cooper-pairs 
which have one spin-up electron on sublatticc a, and one spin-down electron 
on sublattice b. In the partially polarized anti-ferromagnet the amplitude 
of the S=l contribution is different from the S=0 weight a Q = (tt, tt). The 
degree of polarization of the Cooper-pairs depends on the details of the pairing- 
mechanism and the degree of polarization of the underlying anti-ferromagnetic 
state. For vanishing spin-polarization SU(2)-symmetry should be restored, and 
the Cooper-pairs are either in an S=l state or in an S=0 state, depending on 
the parity of Vk/uk- 



5 Supersymmetric rotations and SO(5) 

A different situation arises, when the anti-ferromagnetic correlations and the 
superconducting correlations are stabilized on comparable energy scales. In 
this case the anti-ferromagnetism can be partially or totally suppressed when 
superconducting order sets in and vice-versa. An elegant approach to this 
phenomenon is the group-theory, which uses the super-symmetric S0(5) exten- 
sion of the direct product of U(l) corresponding to the superconducting phase, 
and of SO (3) corresponding to the anti-ferromagnetic order parameter'^^. In 
this so-called SO (5) theory, the anti- ferromagnetic state and the d-wave su- 
perconducting state are regarded as two different projections of a generalized 
higher dimensional order parameter. It is tempting to regard Eq. 13 as a 
manifestation of precisely this SO (5) symmetry, namely, that the d-wave su- 
perconducting condensate can be transformed into anti-ferromagnetic state by 
adiabatically adding the additional triplet condensate at G = (tt.tt). However, 
in Eq. 13 the singlet at G = has s-wave symmetry if the anti-ferromagnetic 
state is of the conventional variety. Hence, although our wave-function does 
support a super-symmetric rotation from anti-ferromagnetism to s-wave pair- 
ing by mixing in a triplet-condensate at G = (tTjTt), a similar construction 
connecting d-wavc pairing to conventional anti-ferromagnetism appears to be 
absent?". In principle we could postulate an unconventional anti-ferromagentic 
order parameter /x^ which has a d-wave symmetry with nodes along the di- 
agonal directions. Such an order parameter transforms to a d-wave supercon- 
ducting order parameter if we reduce the singlet and triplet condensates at 
G = (tt, tt) to zero. 
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Figure 2: Non-BCS scenario for the evolution from the normal to the superconducting state. 
Panels a and b correspond to a situation where pairs axe already formed, but no condensation 
has yet occurred. The superconducting state is only reax;hed when the i5-function gains a 
finite intensity, as indicated in panel c. 



6 Local pairs 



It has often been speculated that a remnant of the Cooper-pairs may exist even 
if the material is not superconducting. This would require that those pairs are 
not condensed in the same state. The variational wave-function proposed by 
Bardeen, Cooper and Schrieffer does not include this possibility. On the other 
hand the situation sketched in Fig. la can be easily generalized to describe 
a paired state with no long range superconducting order, by replacing the 
Dirac (5-function of the center-of-mass coordinate with the Lorentzian 1/(1 + 
l^|Gp). The width of this distribution, l/l then corresponds to the inverse 
of length over which the pairs maintain the phase-coherence. A transition to 
the superconducting state as a function of temperature, or another control 
parameter, is in this scenario characterized by the appearance of the Dirac S- 
function at the origin. Note, that for the existance of superconductivity it is not 
necessary that the entire distribution collapses into the (5-function, not even at 
T = 0. This scenario, depicted in Fig. 2, is reminiscent of the experiments with 
Bose-Einstcin condensation. Yet, the pairs in this wave- function are not real 
bosons. In principle it is not even required that they are strongly localized in 
space, although the conditions giving rise to this type of phase transition would 
typically cause the pairs to be rather small. States of the type displayed in Fig. 
2a have a simple representation in real space. Starting from Eq.2, we insert 
the real-space representation of the creation operators cj[.^ = J2n c'^^^c]^ , and 



8 



obtain 

|*2;v) = C{^firm,rn)cl^cl^r\0) (16) 

m,n 

where f{rm,rn) = J2k G e''^'''^^'"~^"-^~^^'^'^^ f{k, G) is the double Fourier-transform 
of the momentum distribution function. A state with near-neighbor resonating 
bonds is constructed by defining 

firm, Tn) = e-|'^-+'^--«°l/(20<5<„_„> {(x^ - xr,) + i{ym - Vn)}'' (17) 

where I measures the range of the RVB correlations measured from an arbitrary 
reference point (Rq) where the correlations are maximal, (^<r?T,,n> 

selects out 

nearest-neighbour bonds, and the last factor selects the state with relative 

angular momentum L. A Gutzwiller projection is still required if we want to 
exclude the doubly occupied sites. The representation in k-space is 

G.k 



where the +/— sign refers to L = or L = 2 angular momentum states 
respectively, 4'{G) is a random phase, and P is the Gutzwiller projection op- 
erator. The true RVB stat^'*'^'®'^ corresponds to taking the limit Z — > oo in 
the above expression. In this limit this state becomes equivalent to a BCS- 
type superconducting wave-function with either extended s-wave symmetry or 
d-wave symmetry, with the two paired holes on a nearest-neigbour distance. 
In a superconductor with a small gap, the fc-dependent factor inside the curly 
brackets has a strong k-dependence, corresponding to singlet-bonds on a range 
much longer than a nearest neighbor distance in /(rm,r„). 



7 Finite range magnetic correlations 

Likewise a wavefunction corresponding to a state of short-range magnetic cor- 
relations can be easily constructed by replacing the Dirac (5-functions in Fig. 
2b with distribution functions with a finite width. The corresponding struc- 
ture sketched in Fig. 3a corresponds to the magnetic structure-factor mea- 
sured with neutron scattering, and the width of the distribution is inversely 
proportional to the anti-ferromagnetic correlation length. Like for the super- 
conducting state, long-range anti-ferromagnetic order is characterized by the 
condensation of part of the spectral weight into the (5-function at G = (7r,7r). 
By shifting the peak in the distribution function from (tt, tt) to another value of 
the wavevector G, the formalism can be used to describe short- or long-range 
incommensurate spin-density waves. 
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Figure 3: Scenario for the evolution from a state with short range anti- ferromagnetic order 
(panel o)to the Neel state with long range order (panel 6). 



8 Summary and conclusions 

We have generalized earlier expressions of the many-body wavefuctionf'*'^^ to 
a formula, Eq.l, which envelopes many known or suspected types of spin- and 
pairing order of the high Tg cupratcs and of some of the heavy fermion su- 
perconductors. We have used this approach to show that anti-ferromagnetism 
can be regarded as the simultaneous condensation of S=0 singlets at G = 
and G = Q, and of S=l triplets at G = Q where Q is the anti-ferromagnetic 
Bragg- vector. The possibility of superconductivity and anti-ferromagnetism 
occurring simultaneously can be easily described this way. In principle the 
calculation of matrix-elements of the Hamiltonian, as well as most operators 
corresponding to experimentally measurable quantities, is unproblematic at 
least in the examples where the ccntcr-of-mass momentum, G, has a single 
discrete value. The general case presented by the proposed wave-function al- 
lows G to be distributed over a finite range of fc-space. In the latter case it 
appears to be more problematic to obtain analytical expressions of the afore- 
mentioned matrix-elements. The latter forms, as yet, the main obstacle for 
using Eq. 1 in a variational calculation of the ground state and the excita- 
tion spectrum. If this obstacle can be taken, this could offer a straightforward 
analytical approach to the rich phase diagram of the high Tc cuprates. 

The work presented in this paper was developed during a long term col- 
laboration with M. J. Rice. His pedagogical teachings have largely influenced 
the present paper, written in his memory. It has been a great privilege to 
collaborate with Michael. 
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